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Abstract 

Detection of relic gravitational waves, through their imprint in the cosmic microwave background radiation, is one of the most 
important tasks for the planned CMBPol mission. In the simplest viable theoretical models the gravitational wave background 
T-H , is characterized by two parameters, the tensor-to-scalar ratio r and the tensor spectral index nj. In this paper, we analyze the 
T-H ■ potential joint constraints on these two parameters, r and nt, using the potential observations of the CMBPol mission, which 
' is expected to detect the relic gravitational waves if r > 0.001. The influence of the contaminations, including cosmic weak 
' lensing, various foreground emissions, and systematical errors, is discussed. 
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I. INTRODUCTION 



The detection of primordial gravitational waves is rightly considered a highest priority task for the upcoming 
observation missions A stochastic background of the relic gravitational waves (RGWs), produced in the very 
early Universe due to the superadiabatic amplification of zero point quantum fluctuations of the gravitational field, 
is a necessity dictate by general relativity and quantum mechanics Q . So the detection of RGWs maybe provide the 
unique way to study the birth of the Universe, and test the applicability of general relativity and quantum mechanics 
in the very high energy scale [3] . 

In a whole range of scenarios of the early Universe, theprimordial power spectrum of the RGWs can be well described 
by a power-law form in a fairly large frequency range [^-l2|- Thus, the RGW backgrounds are conventionally simply 
characterized by two parameters, the so-called tensor-to-scalar ratio r and the primordial power spectral index of 
RGWs Tit , where r describes the amplitude of the primordial spectrum, and nt denotes the tilt of the spectrum. So 
the constraints on r and rit will give us a direct glimpse into the physical conditions in the early Universe. In particular, 
they will allow us to place the constraints on the Hubble parameter of the early Universe, and time evolution of this 
Hubble parameter. Unfortunately, the models of the early Universe cannot give a definite prediction for the values of 
r and rit, i.e., diff'erent infiationary models predict the quite different values of r and rit llOll . especially some string 
motivated inflationary models predict a very small gravitational waves with r ^ 10"''' So, the only way to 

determine them is by the observations. 

The RGWs leave well understood imprints on the anisotropics in temperature and polarization of cosmic microwave 
background radiation (CMB) jT^ - [T8l |. More specifically, RGWs produce a specific pattern of polarization in the 



CMB known as the _B-mode polarization [13[. Moreover, RGWs produce a negative cross-correlation between the 



temperature and polarization known as the TiJ-correlation at the low multipoles £ < 50 [17|, |18[ . The theoretical 
analysis of these imprints along with the data from CMB experiments allows to place constraints on the parameters 
r and nt describing the RGW background, which provides the unique way to detect the RGWs at the very low 
frequencies (lO-^"^ - IQ-^^ Hz). 

The current CMB experiments are yet to detect a definite signature of RGWs [l^, although a hint of RGWs is 
found in the WMAP data [13] . A number of authors have discussed the possibility of RGW detection by the launched 
Planck satellite [20l - |23| . The results show that, due to the fairly large instrumental noises, only if the tensor-to-scalar 
ratio is large (r > 0.05), the Planck satellite is expected to have a detection. In the previous paper [2J|, we also found 
that the Planck mission cannot give a good constraint for the spectral index nt , even if the tensor-to-scalar ratio is as 
large as r = 1. In addition, various ground-based [1^1 and balloon-borne [IB] CMB experiments are expected to have 
the better detection abilities, which can constrain the parameter r fairly well when r > 0.01. However, an accurate 
constraint of nt is still unexpected, due to the small partial sky observation or the short time observation i27j . 

The accurate determination of the RGWs requires the full sky and the long time observations by the CMB exper- 
iment with the quite small instrumental noises. The future space-based mission focus on the CMB polarization [2^ 
(here and in the following, we use the label 'CMBPol' to refer it) provides an excellent opportunity to realize it (the 
similar projects, such as B-Pol [2^ and LiteBird [3^, are also proposed). The instrumental noises of CMBPol mission 
are more than 100 times smaller than those of Planck mission. If the foreground contaminations and the systematical 
errors can be well controlled, the signature of RGWs can be well detected, as long as r > 0.001 [11]. This will provide 
an observational tool to distinguish the different inflationary-type models. 

In this paper we shall analyze the joint constraints on two parameters r and nt that would be feasible with the 
analysis of the observations from the planned CMBPol mission. We shall detailedly discuss the constraints of r, nt 
and the best-pivot wavenumber kt depending on the input (or true) value of the tensor-to-scalar ratio r. We discuss 
the effects of various contaminations from the cosmic weak lensing, foreground radiations and the beam systematics. 

The outline of the paper is as follows. In Sec. |TT]we shall introduce and explain the notations for the power spectra 
of gravitational waves, density perturbations and various CMB anisotropy fields. Furthermore, in this section, we 
shall briefiy introduce the existence of the best-pivot wavenumber for the detection of RGWs in the CMB. The 
analytical formulae and explanation for the fc^ , Ar and An* will also be discussed. In Sec. IIIIl by using the analytical 
formulae and only considering the instrumental noises, we shall discuss the values of fcj , Ar and Arii for different 
input (or true) value of r. Sec. IIVI is contributed to show the effect of the cosmic lensing contamination, and Sec. 
IVl is contributed to show the effect of the foreground radiations contamination. In Sec. I VII we discuss the effects of 
various beam systematics for the determination of the parameters r and nt- We also discuss the requirement of the 
CMBPol's systematics, if the biases of the parameters r and nt are ignorable. Finally, Sec. I VIII is dedicated to a brief 
discussion and conclusion. 
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II. OPTIMAL PARAMETERS AND THEIR DETERMINATIONS 



The main contribution to the observed temperature and polarization anisotropics of the CMB comes from two types 
of the cosmological perturbations, densi ty p erturbations (also known as the scalar perturbations) and RGWs (also 
known as the tensor perturbations)! a Il2l . IT3I I , which are generally characterized by their primordial power spectra. 
These power spectra are usually assumed to be power-law, which is a generic prediction of a wide range of scenarios of 
the early Universe, including the inflationary models. In general there might be deviations from a power-law, which 
can be parameterized in terms of the running of the spectral index (see for example [lol.[3ll|). but we shall not consider 
this possibility in the current paper. Thus, the power spectra of the perturbation fields have the form 

Ps{k) - AAko){k/kor^-\ Pt{k) = At{ko){k/kor\ (1) 

for density perturbations and RGWs respectively. In the above expression ko is an arbitrarily chosen pivot wavenum- 
ber, Us is the primordial power spectral index for density perturbations, and rit is the primordial power spectral index 
for RGWs. As{ko) and At{ko) are normalization coefficients determining the absolute values of the primordial power 
spectra at the pivot wavenumber fco . The choices of = 1 and rit — correspond to the scale- invariant power spectra 
for density perturbations and gravitational waves respectively. 

The relative contribution of density perturbations and gravitational waves is described by the so-called tensor-to- 
scalar ratio r, which is defined as follows 

As[ko) 

Note that, in defining the tensor-to-scalar ratio r, we have not used any inflationary formulae which relate r with 
the physical conditions during inflation and the slow- roll parameters (see for example [12]). Thus, our definition 
depends only on the power spectral amplitudes of density perturbations and RGWs, and does not assume a particular 
generating mechanism for these cosmological perturbations. The RGW amplitude At (/cq) = ''(^0)^5(^0) provides us 
with direct information on the Hubble parameter of the very early Universe [31| . More specifically, this amplitude is 
directly related to the value of the Hubble parameter _ff at a time when wavelengths corresponding to the wavenumber 
fco crossed the horizon 33.1 



(3) 

ko I a—H 



where Mpi = I/VSttG is the reduced Planck mass. If we adopt As = 2.445 x 10 ^ as predicted by the WMAP5 
observations [s^l, the Hubble parameter is if ~ 2.67r^/^ x lO^^'^GeV, only depending on the value of r. In the canonical 
single- field slow-roll inflationary models, the Hubble parameter directly relates to the energy scale of inflation y^^^. 
The relation ([3]) follows that V^^'^ ~ 3.35r^/* x lO^^GeV, which has been emphasized by a number of authors. 

Assuming that the amplitude of density perturbations ^^(fco) is known, taking into account the definitions ([1]) 
and ([2]), the power spectrum of the RGW field may be completely characterized by tensor-to-scalar ratio r and the 
spectral index nt. It is important to mention that, for spectral indices different from the scale-invariant case (i.e., 
when ris 7^ 1 or/and rit 7^ 0), the definition of the tensor-to-scalar ratio depends on the pivot wavenumber fcp. If we 
adopt a different pivot wavenumber fci, the tensor-to-scalar ratio at this new pivot wavenumber r(fci) is related to 
original ratio r(fco) through the following relation (which follows from the definitions ([T]) and ([2|)) 

fcl 



r(fci) = r(fco)(^^j . (4) 

Let us now turn our attention to CMB. Density perturbations and gravitational waves produce temperature and 
polarization anisotropics in the CMB, which are characterized by four angular power spectra Cj, , Cf and Cf as 
functions of the multipole number I. Here Cj is the power spectrum of the temperature anisotropics, Cf and Cf are 
the power spectra of the so-called E and B modes of polarization (note that, density perturbation do not generate 
_B-mode of polarization !13i), and Cf is the power spectrum of the temperature-polarization cross correlation. 

In general, the power spectra Cj (where Y = T,E,B or C) can be presented in the following form 

= Cf(dp)-fCf(gw), (5) 
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where C^(dp) is the power spectrum due to the density perturbations, and C^(gw) is the power spectrum due to 
RGWs. In the case of RGWs, the CMB power spectra can be presented in the following form [l^. [isj 



2 

, for y = T, E, B, 

Cf (gw) = (AnfJ fp.ik) [A(.7(fc)Aj;(fc)] . (6) 

The transfer functions A^'' (fc) (see [iBl for details) in the above expressions translate the power in the metric 
fluctuations (gravitational waves) into corresponding CMB power spectrum at an angular scale characterized by 
multipole £. In this work, for numerical evaluation of the CMB power spectra due to density perturbations and 
gravitational waves, we use the publicly available CAMB code [35| . 

Since we are primarily interested in the parameters of the RGW field, in the analytical and numerical analysis 
below we shall work with a fixed cosmological background model. More specifically, we shall work in the framework 
of ACDM model, and keep the background cosmological parameters fixed at the values determined by a typical model 

M 

h ^ 0.705, nbh^ = 0.02267, f^^/i^ = 0.1131, = 0, Tre^on = 0.084, A, = 2.445 x 10"^ (7) 

Furthermore, the spectral indices of density perturbations and gravitational waves are adopted as follows for the 
simplicity. 

Us = 1, nt = 0. (8) 

Note that throughout this paper, we have considered the simplest cosmological model. In the more general con- 
sideration, one should also include the running of the spectral indices [Slj , the details of the reionization history [s^ 
and so on, which have been ignored in this paper. 

The CMB power spectra Cj are theoretical constructions determined by ensemble averages over all possible real- 
izations of the underlying random process. However, in real CMB observations, we only have access to a single sky, 
and hence to a single realization. In order to obtain information on the power spectra from a single realization, it 
is required to construct estimators of power spectra. In order to differentiate the estimators from the actual power 
spectra, we shall use the notation Dj to denote the estimators while retaining the notation Cj to denote the power 
spectrum. It is important to keep in mind that the estimators Dj are constructed from observational data, while the 
power spectra Cj are theoretically predicted quantities. The probability distribution functions for the estimators are 
described in detail in [18| (see also ,2^ iST,, ,38] ) , which predicts the expectation values of the estimators 



i^n-Q, (9) 



and the standard deviations 



^""^^^ = (2£+l)/3., ' 

where /sky is the sky-cut factor. In this paper, we use /sky — 0.8 for the CMBPol survey. are the noise power 
spectra, which are all determined by the specific experiments. In this formulas, the possible bias generated by the 
beam systematics has not been considered (see Sec. IVII for details). 

In order to estimate the parameters r and Ut characterizing the RGW background, we shall use an analysis based on 
the likelihood function [39]. The likelihood function is just the probability density function of the observational data 
considered as a function of the unknown parameters (which are r and rit in our case). Up to a constant, independent 
of its arguments, the likelihood function is given by 

C^Y[f{Df,Dj,Df,Df), 

I 

where the function f{D9, Dj , D? , D?) is explained in detail in the previous works p2l. [23|. 

In the previous wor ks f23 .l24|. we have discussed how to constrain the parameters of the RGWs, r and nj, by the 
CMB observation. In 24|, we found that in general, the constraints on r and rit correlate with each other. However, 
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if we consider the tensor-to-scalar ratio at the best-pivot wavenumber k^, i.e. r* = r{k^), the constraints on r and rit 
becomes independent of each other, and the uncertainties Ar and Ant have the minimum values. We have derived the 
formulae to calculate the quantities: the best-pivot wavenumber k^, and the uncertainties of the parameters Ar and 
Ant, which provides a simple and quick method to investigate the detection abilities of the future CMB observations. 
We shall briefly introduce these results in this section. 
It is convenient to define the quantities as below, 



where Cj{g-w) is the CMB power spectrum generated by RGWs, and cr^i- is the standard deviation of the estimator 

DJ , which can be calculated by Eq. ([TU| . We should notice that, the quantity c?^' is dependent of random date 
DJ . By considering the relations in © and ([S]), we can obtain that {d^) = a|^, which shows that <f( is an unbiased 
estimator of aj . it is the so-called best-pivot multipole, which is determined by solving the following equation [2^]: 



EE' 



arK = 0. (12) 

So the value of depends on the cosmological model, the amplitude of RGWs, and noise power spectra by the 
quantity aj . The best-pivot wavenumber kt relates to £t by the approximation relation p^ . 

k; ~ e*t X 10~*Mpc-^ . (13) 

The numerical factor here mainly reflects the angular-diameter distance to the last scattering surface. 

Once the value of it is obtained, the uncertainties Ar* and Ant can be calculated by the following simple formulae 



- jEE«r' - V jEE(«r^.*)'- (i^) 



^ I Y ^ e Y 

As usual, we can define the signal-to-noise ratio S/N = r* / Ar* . Using p4)) . we get 



^/^ = wEE«r- (15) 

\ I Y 

In the previous work j24j . we found the uncertainty of r{kQ), the tensor-to-scalar ratio at the pivot wavenumber 
fcg kt , is larger than Ar*. The value of Alnr(fco) is fairly well approximated by the following formula 



Alnr(fco) = y/{Ar*/r*)^ + (ln(fco/fc*)Ant)2. (16) 

The smallest uncertainty on tensor-to-scalar ratio r is achieved for the choice of the pivot scale at feg = fc^*. This 
justified the title 'best' pivot wavenumber for k* . We should notice that the values of k't, S/N, Alnr and Ant only 
depend on the input (or true) cosmological model, but not on the data Dj. In Fig. [H we plot the value of A Inr(fco) as 
a function of the pivot scale ka, where the input model has r — 0.3, and the Planck instrumental noises are considered 
(see [l^l for details). As expected, when fco ^ or ko <^kl, the uncertainty becomes much larger than Ar*. 

The likelihood function in ([TT|) has the maximum value at (rJ^L, n^ML)- The values of r^^ ^^id njML depend on the 
data DJ , different from the value of Ar* and Ant- In the previous work [23|, we found that the values of r^L ^''^d 
JT-tML can be very well approximated by the follows 

* * Til YIy ^J'^J ^ I HiYIy '^J '^J^l /lyN 

ryih-r ^ ^ ^y2 ' '^tML - + („Yi,*-.2 ' (1") 

which depend on the data by the quantity d^. If the CMB estimator Dj is unbiased for Cj , as discussed above, we 
have (r^j^) = r* and {ntuh) = n-t, where Eq. is used. These show that r^^ ^-nd ntuh are the unbiased estimators 
for r* and rif, respectively. However, when Dj is a biased estimator for , ''ml ^'^d rif ml will also be the biased 
estimators for r* and nt, respectively (see Sec. IVII for details), which brings the errors for the detection of RGWs. 

The detection ability of the CMB experiment strongly depends on the noise levels, which include the instrumental 
noises, cosmic lensing contaminations, foreground radiation contaminations and the beam systematics. In the following 
sections, we shall discuss these effects separately. In addition, due to the partial sky survey, the leakage from the 
i?-polarization into the B-polarization could be another kind of contamination. However, it was found that, this E-B 
mixture can be properly avoided (or deeply reduced) by constructing the pure i?-mode and i?-mode polarization fields 
(see [40l - |4^ for details). So we shall not discuss this topic in this paper. 
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FIG. 1: The uncertainty of r(fco) for the different pivot wavenumber ko- 



TABLE I: Experimental specifications for the mid-cost (EPIC-2m) CMBPol mission and the low-cost (EPIC-LC) CMBPol 
mission. [281 





Frequency [GHz] 


45 


70 


100 


150 


220 


EPlC-2m 


[arcmin] 


17 


11 


8 


5 


3.5 




At [/^K-arcmin] 


5.85 


2.96 


2.29 


2.21 


3.39 




Frequency [GHz] 


40 


60 


90 


135 


200 


EPIC-LC 


[arcmin] 


116 


77 


52 


34 


23 




At [/^K-arcmin] 


15.27 


8.23 


3.56 


3.31 


3.48 



III. CMBPOL INSTRUMENTAL NOISES' CONTAMINATIONS 



In this section, we shall discuss the determination of RGWs, when only taking into account the instrumental noises 
of the CMBPol mission. 

For a single frequency channel i, we assume Gaussian beams. The noise power spectrum (after deconvolution of 
the beam window function) is 



iVi^,,« = (AT)2exp 



8 In 2 



and 



8 In 2 



(18) 



(19) 



where Op is the full width at half maximum (FWHM) of the beam i. At and Ap are the noises for the temperature 
and polarizations, which relate by Ap = \/2At. The values of the A^ and Ap depend on the number of the detectors, 
the integration time and the survey area. 
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FIG. 2: The instrumental noise power spectra N^^ ^ of EPIC-2m (black line) and EPIC-LC (red line). 

If the experiment inchides several different channels, we need to generalize the above considerations. The optimal 
channel combination of these channels gives the total effective instrumental noise psl] , 

[NL,ir' = Y.[^Ls)V'^ (20) 

i 

where i runs though the channels, Ny^f{i) is the instrumental noise bias of the channel i^^. In this section, we shall 
only consider the instrumental noises, i.e. 

Nl ^ A^i^M- (21) 

Since the precise experimental specifications of CMBPol have not yet been defined, we will consider two different 
cases (EPIC-2rii) and (EPIC-LC) suggested by [11] (recently, an EPIC-Intemediate Mission is also suggested by the 
CMBPol team (ijl)- The experimental specifications are given in Table HI where 2-year design life is assumed [61j. In 
Fig. [21 we plot the polarization noise spectra N^^ ^ of EPIC-2m and EPIC-LC, respectively For EPIC-2m, when 
£ < 100, we find that N^f^^ ^ ^ 2.7 x 10~^/iK^, which is nearly 400 times smaller than that of the Planck mission (7 
frequency channels from SOGHz to 353GIIz and 28-month surveying time are assumed "21]). Even for the EPIC-LC, 
when I < 100, we have ^ 6.2 x 10^^/iK^, 200 times smaller than that of the Planck mission. So comparing 

with Planck mission, CMBPol is much more sensitive for the detection of CMB polarization. From Fig. [H we also 
find that even for the model with quite small r = 0.01, the value of N^f^^ ^ is smaller than that of Cf when I < 120 
for EPIC-2m, and when £ < 80 for EPIC-LC. So the CMBPol mission provides an excellent opportunity to detailedly 
observe the peak of at € ~ 80. 

Let us discuss the constraints on the gravitational waves by the potential observations of CMBPol mission. We 
shall discuss the values of the best-pivot scale fcj , the signal-to-noise ratio S/N and the uncertainty of the spectral 
index Ant, by considering the CMBPol instrumental noises. 

The value of fc^* directly relates to the best-pivot multipole £j by Eq. (fT3)) . and the value of is obtained by solving 
the equation in (IT^ . By using (l?T|) . we obtain the value of as a function of the input (or true) value of the 
tensor-to-scalar ratio r* for EPIC-2m and EPIC-LC, which are plotted in Fig. [3] (left panel). We find that, in both 
cases, the value of becomes larger with the increasing of r* . For EPIC-2m, we have = 43 for r* = 0.001, and 



7 




FIG. 3: The figures show the values of the best-pivot multipole i* (left panel), signal-to-noise ratio S/N (middle panel) and 
the uncertainty of the RGW spectral index Ant (right panel) as functions of r*. The black solid (dash-dotted, dashed) lines 
correspond to the EPIC-2m instrumental noises (instrumental noises -I- the reduced cosmic lensing contaminations, instrumental 
noises -I- cosmic lensing contaminations), and the red solid (dashed) lines correspond to the EPIC-LC instrumental noises 
(instrumental noises -|- the cosmic lensing contaminations). 



£* = 137 for r* — 0.1. For EPIC-LC, the value of is smaller than that of EPIC-2in, due to the larger noise level 
and the larger beam FWHM of EPIC-LC. When r = 0.001, we have £; = 26, and when r = 0.1, we have 4* = 87. 
These reflect that gravitational waves in the frequency range k ~ O.OlMpc"^ will be best constrained by the future 
CMBPol observations, unless the value of r is extremely small. This is because the main contribution comes from the 
observation of the peak of B-polarization at £ ^ 80. We should remember that this is different from the Planck case, 
where ^ 10, due to the main contribution of the reionization peak of B-polarization [isl [2^ [2^. 

The signal-to- noise ratio is calculated by Eq. ([15]). By using Eq. (|2ip . we get the value of S/N as a function of 
r* for both EPIC-2m and EPIC-LC, which are shown in Fig. [3] (middle panel). As expected, the signal of RGWs 
can be very well determined by the CMBPol mission. Even for the model with r = 0.001, we can have S/N — 19 for 
EPIC-2m and S/N = 10 for EPIC-LC, when only considering the corresponding instrumental noises. When r ~ 0.1, 
we have S/N = 122 for EPIC-2m and S/N = 81 for EPIC-LC. 

We can also calculate the value of Ant by using Eq. ([T4)) and the value of £1 given in left panel of Fig. [3] The results 
are shown in Fig. [3] (right panel). As expected, the value of Ant decreases with the increasing of r* . For EPIC-2m, 
we have An* = 0.06 for the model with r* = 0.001, and Ant = 0.01 for the model with r* — 0.1. This uncertainty is 
about 20 times smaller than that given by Planck satellite [l^, H^l- This constraint, combining with Ar* , will give a 
quite sensitive way to differentiate various inflationary type models. For the EPIC-LC, the uncertainty of nt is about 
2 times larger that of EPIC-2m. When r = 0.001, we have Ant — 0.10, and when r = 0.1, we have Ant = 0.02. 

It is necessary to discuss the contributions of S/N and Ant from every multipole, which can be very easily analyzed 
by the analytical formulae. From Eqs. (IT4l) and (fT5|) . we find these two quantities can be rewritten as follows 

(W-EE^r- (l/An,)^=EE(«F^^*)'- (22) 

e Y e Y 

They are the simple sums of the contributions from each multipole £ and CMB information channel Y . We plot the 
functions of J^y ^^'^ J^Yi'^JK)'^ ^ function of £ for two different models (r = 0.01 and r = 0.1). The results 
are shown in Fig. 21 Left panel shows that all these four lines are peaked at £ ~ 100, which is close to the peak 
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FIG. 4: The figures show the values of af^ (left panel) and '^yi'^^ ^e)^ (right panel) as functions of multipole i for the cases 
with the different r*. Ifere, we have not considered the contaminations from cosmic weak lensing and foreground emissions. 



of i?-polarization. This reflects that when r > 0.01, the main contribution comes from the observation in the range 
i ~ 100, consistent with our previous discussion. This is different from the case of Planck satellite [H, Hljj where the 
reionization peak at ^ ~ 6 is extremely important. In the CMBPol case, the contribution from the largest scale £ < 20 
is unimportant due to the cosmic variance, and the contribution from the small scale £ > 300 is also unimportant for 
the large instrumental noises. However, it is important to mention that if r ^ 0.01, similar to Planck satellite, the 
reionization peak at ^ < 10 again becomes the main contribution for the total S/N . 

However, it is different for ^y('^JK)'^' which stands for the individual contribution for l/An^. From the right panel 
of Fig. m we find that this function is sharply peaked at the largest scale £ < 30, and the contribution from intermedial 
scale around the best-pivot multpole is very small. These can be easily understood, the quantity bg = lii(£/£l) is zero 
when £ = £f, which follows that J^yi^I^l)'^ — a,t £ — £1- Only if £ ^ or ^ > £1, b} has a large value, and follows 
a large X^yC^^^I)^- Especially, the contribution from £ <^ is very important. For example, when £^ = 137, 61^2 
30 times large than bg^^^Q. This reflects that the constraint on the tilt of the primordial gravitational waves power 
spectrum strongly depends on the observations in a large scale range. The cosmic reionization is very important for 
the constraint of fit, although it might not be so important for the constraint of r for the CMBPol observations. 



IV. COSMIC WEAK LENSING CONTAMINATION 

In [l^, it was pointed out that the gravitational waves result in CMB polarization with a B-mode, whereas density 
perturbations do not. Thus, the signal of gravitational waves could not be confused with density perturbations by 
detecting the S-polarization. Although, the amplitude of the B-polarization is expected to be quite small, it gives a 
clear information for gravitational waves. However, when taking into account the second-order effect, the i?-mode can 
also arise from the lensing of the E'-mode by density perturbations along lines-of-sight between the observer and the 
last-scattering surface [44|. The scalar contribution to the i?-mode power spectrum is shown in Fig. [5] (grey dashed 
line). When £ < 200, it is nearly a white spectrum with the amplitude C^^^ ^ ~ 2 x 10~^^K^, which is 7 times larger 
than the instrumental noises of EPIC-2m, and 3 times larger than that of EPIC-LC. 

When the instrumental noise of the CMB experiment is sufficiently small, as the CMBPol mission, the gravitational 
lensing contribution to the large-scale _B-mode becomes one of the limiting sources of contamination for constraining 
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FIG. 5: For EPIC-2m, the figures show the values of (l^ft panel) and X]y ('^F^D'^ (right panel) as functions of multipole 

i for the cases with and without the reduced cosmic lensing contamination with the residual factor a'°"'' = 0.5. 

the RGWs. High-sensitivity measurements of smah-scale _B-modes can reduce this contamination through a lens 
reconstruction technique, which has been discussed by a number of authors (see for instance [45l - l47j ). 

The effect of cosmic lensing contamination for the detection of RGWs can be easily discussed. The reduced lensed 
B-mode polarization can be treated simply as a well-known noise for gravitational waves in the likelihood analysis, 
i.e. 

^A^ifM + Qfns,.xa'°-, (23) 

where we have defined the residual factor a^^™ for the lensed _B-polarization. Note that in the real situation, the 
lens-induced i?-modes are non-Gaussian, so we should not behave exactly it as the additional Gaussian noise (see [i^l 
and references therein). However, on the scales relevant for B-mode detection in this paper, the non-Gaussianity has 
only a minor effect, which has been ignored in our discussion. In general, we have a^''^^ < 1^ with the equality holding 
for the lensed B-mode is not reduced. The reduction of gravitational lensing contribution strongly depends on the 
instrumental noise level, the beam FWHM, the foregrounds and the instrumental systematics. In the work [48], the 
authors found that, based on the noise level of EPIC-2m, one can expect to have cr''^"^ ^ 0.5. However, for EPIC-LC, 
it is very difhcult to reduce the cosmic lensing contamination due to the large beam FWHM. We should mention that 
since the value of C^^^^ ^ is much larger than the instrumental noises of CMBPol mission, in the total effective noise 
Nf ^ the cosmic lensing contamination becomes the dominant portion. 

We have calculated the constraints of the gravitational waves by taking into account the cosmic lensing contam- 
inations. The values of ^J" , SjN and Ant are shown in Fig. [3l where a^™^ = 0.5 and a^^^^ = 1 are considered for 
EPIC-2m, and cr''^"*' = 1 is considered for EPIC-LC. The left panel shows that, the best-pivot multipole is shifted to 
smaller scale by the lensing contamination. We find the value of SjN is much reduced by the lensing contamination, 
especially for the case with small tensor-to-scalar ratio. When r — 0.001, we have SjN = 7 for EPIC-2m (with 
^lons _ Q g-j g^j-^jj gjj^ _ ^ epjc_lc^ which are much smaller than the corresponding values with only instrumental 
noises. The uncertainty of nt is also much increased by the cosmic lensing, especially for the case with small r. When 
r = 0.001, EPIC-2m with cr^™" = 0.5 can give An* = 0.12, and EPIC-2m with cr^™'' = 1 can give An* = 0.16, which 
are more than 2 times larger than those in the case without cosmic lensing, and are fairly loose to differentiate various 
inflationary models. However, when the tensor-to-scalar ratio is r = 0.1, we have Arii = 0.02 for EPIC-2m, which is 
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TABLE II: Assumptions about foreground emissions. [2^ 
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FIG. 6: The total effective noise power spectra N^g , 
EPIC-2m and the right panel is for EPIC-LC. 



when considering the foreground contaminations. The left panel is for 



still a quite tight constraint. 

We have also investigated the contribution of every multipole for S/N and An*. The results can be found in Fig. [SJ 
"where we have focused on the EPIC-2m mission and cr'''"'^ = 0.5 is used for the case with cosmic lensing contamination. 
We find that the peak in each case is much reduced by the lensing contamination. 

It is interesting to mention that for the high-sensitivity detectors the residual lensing noise dominates over the 
instrumental noises, and place the detection limit for CMB experiments [4!||. In ^3], the authors claimed that, for a 
extreme high-sensitivity detector, a reduction in lensing po"wer by a factor 40 is possible using approximate iterative 
maximum-likelihood method. If we consider this residual as the lo"wer limit of the reduced lensing noises, we find that 
r > 3.7 X 10^^ can be detected at more than 2-cr level in absence of sky cuts, foregrounds and instrumental systematics 
[2^ . This can be treated as the detection limit of the CMB experiments for gravitational waves. This lower limit 
corresponds to the Hubble parameter H ~ 3.1 x lO^^GeV, and the energy scale of inflation T/^/"* ~ 1.5 x lO^^GeV. In 
this limit case, the uncertainty of spectral index also becomes very small. When r — 0.1, we have Ant ~ 0.007 (see 
Fig. 2 in [23| for details), placing a very tight constraint on the inflationary models. 



FOREGROUND CONTAMINATIONS 



In this section, besides the instrumental noises and cosmic lensing contaminations, we shall take into account the 
impact of polarized foregrounds on the future CMBPol mission. In this paper, we shall neglect the effect of foregrounds 
on the CMB temperature, as the foreground cleaning is expected to leave a negligible contribution in the temperature 
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FIG. 7: This figure shows the values of £*, S/N and Ant depend on the foreground contaminations for EPIC-2m. 
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FIG. 8: This figure shows the values of £*, S/N and Ant depend on the foreground contaminations for EPIC-LC. 
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FIG. 9: The figures show the values of oj'^ (l^ft panel) and X^vC'^F^^)^ (right panel) as functions of multipole d for the 
cases with different cr*^. 



[50| . CMB polarized foregrounds arise due to free- free, synchrotron, and dust emission, as well as due to extra-galactic 
sources such as radio sources and dusty galaxies. In this paper, we shall only consider only synchrotron and dust 
emission, which are expected to be dominant in the CMBPol frequency range [Slj . 

The synchrotron emission results from the acceleration of cosmic-ray electrons in the ma gnet ic field of Galaxy, 
which has been well measured on large angular scale at 23 GHz by WMAP. Following [2^ [5l|, 52 1, for the frequency 
ly, the scale-dependence of the synchrotron signal may be parameterized as 



Ci/^) = As[- 



2as 



(24) 



The parameters in this formula for the various power spectra are all listed in Table |TT1 where as = —3 is assumed, 
/3^, and (3'^ are the corresponding /3|' for synchrotron emissions. This choice matches the synchrotron emission 
at 23GHz observed and parameterized by WMAP [HS], and agrees with the DASI measurements [s^ . 

Galactic emission in the 100—6000 GHz frequency range is dominated by the thermal emission from warm interstellar 
dust grains. Our knowledge of polarized dust emission is relatively poor, which is expected to be characterized by the 
Planck satellite in the near future. In this paper, we shall adopt the parameterized formula for the dust emission at 
frequency i' as follows, as suggested by [28.. ,52i]. 



2ao 



^hua/kT , 



ohv/kT 



(25) 



where p = 5% and T = 18K. We Ust the other parameters for the various power spectra in Table ITTl 

Various methods have been discussed to subtract the foregrounds by their frequency-dependence (see for instance 
[ssj). In this paper, we shall not discuss the subtraction of the foreground from the signal. Instead, as the previous 
works [H, [111 we assume that the foreground substraction can be done correctly down to a given level, and treat 
these residual foregrounds as a kind of known gaussian noises in the data analysis. However, here we should mention 
that in this case we have assumed one can model and subtract the power spectra of residuals perfectly, and avoid the 
possible issues of bias. This might be a huge challenge for the future polarization observation. 
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If we consider the CMB experiment, including several frequency channels, and the different channels have different 
noise levels, the optimal channel combination gives the effective noise power spectra [2^ [s^l 

= E [KA') + ~' ' (26) 

i 

where i runs though the channels. Ny^g{i) is the instrumental noise power spectra of channel I'i. N^^{i) is the 
residual foreground noises of channel Vi, which is 

NlS)^ E Cl..A^.y^+KA^^)- (27) 

forc=S,D 

Here, C}^^.^ ^ is the model for the power spectrum of the synchrotron and dust signals at the frequency Vi given by Eqs. 
dMl) and (pS)) . and a^^ is the assumed residual factor. Mj^iivi) is the noise power spectrum of the foreground template 
map, as foreground templates are created by effectively taking map differences and thus are somewhat affected by the 
instrumental noise. This term can be calculated by [23, Is^] 

>^lAv^-^ = 7-^ TT77 — ' (28) 



^chan ("chan - l)/4 V^'rcf 

where richan is the total number of channels used, and the reference channel i^icf is the highest and lowest frequency 
channel included in the cosmological analysis for dust and synchrotron respectively, i.e., that listed in Table ID The 
parameters a for the foregrounds under consideration are defined in Table i.e. a = — 3 for the synchrotron 
emissions and a = 2.2 for the dust emission. The quantity N'\^j^ liy-cci) is the white instrumental noise (without the 
beam window function) of the corresponding template channel 52] . 

Thus the total noise power spectra, by combining the multipole-frequency instrumental noises and the residual 
foregrounds, as well as the residual cosmic lensing contamination, are given by 

^ A^cff,£ {X = T, C, E), Nf ^ NSf,i + Q?„s,^ X a'°-, (29) 

where cr'°"'' is the residual factor for cosmic lensing contamination. In this and the following sections, we adopt 
^lens _ Q g EPIC-2m, and cr'°"'' = 1 for EPIC-LC. The effective noise power spectra iV^ ^ strongly depend on the 
residual factor for the foregrounds. When no foreground subtraction is assumed, we have a^^ = 1. In this paper, 
we also consider two assumed residual cases, suggested by CMBPol team [2^: cr^s — g.Ol for the optimistic case, and 
CT^8 = 0.1 for the pessimistic case. In Fig. [SI we plot the effective noise power spectrum ^ with different a^^ for 
EPIC-2m (left panel) and EPIC-LC (right panel). We find that for both EPIC-2m and EPIC-LC, the foregrounds 
increase the effective noise power spectrum in all the multipole range when cr^s = i However, when the foregrounds 
can be well subtracted, the residual foregrounds only increase the noise in the large scale. For cr^s = 0.1, the effective 
noise is increased in the range £ < 300, and for a^^ — 0.01, the effective noise is increased in the range £ < 100. We 
find that even if the optimistic case with a^^ — 0.01 is realized, the effective noise is much larger in the reionization 
peak (i < 20) comparing with the no foreground case. Especially, when r is small, this increased noise is larger than 
the signal Cf, and decreases the contribution of the reionizaiton peak. We have emphasized above, the reionization 
peak is very important for the constraint of spectral index nt for the CMBPol mission, it is predictable that the value 
of Ant would become much larger due to the foreground contaminations, even if the optimistic case is considered. 
This will be clearly shown in the following discussion. 

By using the total effective noise power spectra in ([29]) . we can calculate the values the best-pivot multipole the 
signal-to-noise ratio S/N, and the uncertainty of the spectral index Ant. In Fig. [7]and Fig. [SJ we show the results for 
EPIC-2m and EPIC-LC, respectively. We find that so long as r < 0.3, the value of is increased with the increasing 
of the residual foregrounds. This is because that, the contaminations from foreground are mainly in the low multipole 
(see Fig. [B]). Increasing the foregrounds, the contribution for the detection of RGWs in the low multipoles becomes 
less and less important, induces an increasing of i^. 

From Fig. [7] and Fig. [51 we find that when a^^ = 0.01, the optimistic case we considered, the foregrounds decrease 
the S/N and increase Arij, when the tensor-to-scalar ratio r < 0.03. However, when r > 0.03, the effect of this 
residual foregrounds is negligible. This is also easily understood. The residual foregrounds with cr^s = q.OI only 
increase the total noise power spectra in the largest scale £ < 100. This increased total noises are beyond the signals 
when r is small. We also find that, when cr^s = q.I or a^^ — 1, i.e. the foregrounds are not well subtracted, the effect 
of the foreground contaminations are quite important, especially for the determination of nt- With the decreasing of 
r, the contaminations become more and more important. For the EPIC-2m and the input model with r = 0.1, when 
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a^B = 0.01, we have Ant = 0.020. When a^^ = o.l, the value becomes Ant = 0.024, and when a^^ = i the value 
becomes Ant = 0.045, two times larger than that in the optimistic case. We can investigate another case, for the 
EPIC-2m and the input model with r = 0.01, when a^^ — 0.01, we have Ant — 0.054. When cr^s — q the value 
becomes Ant = 0.098, and when a^^ = the value becomes Ant — 0.330, six times larger than that in the optimistic 
case. So we conclude that if the value of r is not too small, such as r = 0.1, we do not need to remove the foreground 
to a very high level. The difference between optimistic case and pessimistic is very small. However, if the value of r 
is smaller than 0.01, very detailed removal for the foregrounds is very important for the determination of RGWs. 

Fig. [7] and Fig. |5]also show that in the optimistic case, EPIC-2m can detect the signal of RGWs with r — 0.001 
at 5-(T level, and EPIC-LC can detect it at S-cr level. However, in the pessimistic case, EPIC-2m can only detect this 
signal at 2.8-a level, and EPIC-LC can detect it at 1.6-a level. 

As in the previous sections, we can discuss the contribution to the S/N and Ant from the individual multipole, by 
investigating the functions J^y '^J^ ^^"^ l^yi'^J^i)'^- They are plotted in Fig. ^ where we have considered the EPIC- 
2m, and the model with r = 0.1. We find that the foreground contamination mainly affects the S/N by decreasing the 
value of J2y ^J"^ around the peak at £ ^ 100. However, it affects the value of An* mainly by decreasing of J2Yi^J K)"^ 
at the largest scale ^ < 50 and the intermedial range £ ^ 200. 

Now, let us investigate the possible application of the CMBPol mission to differentiate the different inflationary 
models, which plays a role for the future inflation researches. As well known, one of the most important ways to 
distinguish different classes of inflations is to test the so-called inflationary consistency relations. This testing strongly 
depends on the determination of the parameters specifying the relic gravitational waves, i.e. the tensor-to-scalar ratio 
r and the spectral index n*. Now, let us focus on the possible testing of the consistency relation for the canonical 
single- field slow-roll inflationary models [loj . This testing might provide the unique model-independent criteria to 
confirm or rule out this class of models. The possible testing for other inflationary models by the CMBPol mission 
and the ideal CMB experiment can be found in the recent work [56]. 

The consistency relation for the canonical single-field slow-roll infiationary models can be written as [l3l 

r = -8nt. (30) 

We find that this relation only depends on the parameters r and nt. Since the absolute value of nt is expected to 
be one order smaller than that of r, and also the measurement of nt is much more difficult than r, how well we can 
measure the spectral index nt plays a crucial role for testing the consistency relation in Eq. (|30p . 

To access whether CMBPol mission might achieve the consistency relation test goal, in Fig. [7] and Fig. [8] (right 
panels), we compare the values of \nt\ = r/8 with Ant. If Ant < then the constraint on nt is tight enough to 
allow for the testing. From Fig. [71 we find that for the EPIC-2m mission. Ant < \nt\ is satisfied only if r > 0.14 for 
the optimistic case with a^^ = 0.01. In the pessimistic case with a^^ = q jt becomes r > 0.15. Similar results for the 
EPIC-LC can be found in Fig. [H Ant < \nt \ is satisfied only if r > 0.18 for the optimistic case, and r > 0.20 for the 
pessimistic case. So we conclude that the testing of the consistency relation for the canonical single-field slow-roll by 
the CMBPol mission is quite hard. The testing is possible only for some large-field inflationary models. However, we 
should mention that the situation could become quite promising for the general Lorentz- invariant single-field inflations 
and the two-fleld inflations (see [56i] for the details). 

VI. SYSTEMATICS CONTAMINATIONS 

Beyond raw sensitivity requirements for the instrumentals, and the removal of the astrophysical foregrounds, much 
attention has already been given in the literature to the instrumental systematics for the constraints of the cosmological 
parameters and the cosmic weak lensing reconstruction [57l - l60| . The main goal of this section is to illustrate the effect 
of the instrumental systematics and systematically study the impact on the gravitational waves detection for the 
CMBPol mission. 

All the effects of the beam systematics are associated with beam imperfections or beam mismatch in dual beam 
experiments. Several of these effects (e.g. differential gain, differential beam width and the flrst order pointing 
error) are reducible with an ideal scanning strategy and otherwise can be cleaned from the data set. Other spurious 
polarization signals, such as those due to differential ellipticity of the beam, second order pointing errors and the 
differential rotation, persist even in the case of ideal scanning strategy and perfectly mimic CMB polarization. 

The beam systematics due to optical imperfections are dependent of the underlying sky, the properties of the 
polarimeter and the scanning strategy. If the outputs of two beams with orthogonal polarization-sensitive directions 
are slightly different, the temperature anisotropy can leak to the polarization or the i?-mode polarization can leak 
to the i?-mode and vice verse, (see ^5^] for the details). For example, if two beams are exactly same but the overall 
response, this difference of the measured intensity can generate a non-vanishing polarization signal. Another typical 
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example is effect of the beam rotation, which is caused by the uncertainty in the overall beam orientation. This effect 
mixes the Stokes parameters Q and C/, and induces a E-mode and i?-mode leakage. 

The CMB power spectra C^g_^ generated by these systematics are discussed in details by a number of authors. In the 
work [5^ , the authors discussed these effects separately, and got the simple analytical formulas to calculate the leading 
order of the generated power spectra, which are listed in Table IIIIl The formulae in Table IIIII separately describe the 
effects of the following instrumental systematics for experiments with the elliptical gaussian beams: differential gain 
effect, monopole effect, differential pointing effect, quadrupole effect, differential rotation effect. Differential gain can 
induce spurious polarization singles from temperature leakage due to beam mismatch. This effect is described by the 
parameter g = gi — g2, where gi and 52 refer to the gain factors of first and second beams. The differential rotation 
effect is due to uncertainty in the overall beam orientation. This mixes the Q and U Stokes parameters and as a 
result leaks E to B and vice verse. We describe this effect by the parameter e = {ei + S2)/2, where £1 and £2 are 
rotation errors of first and second beams. We note that these two parameters g and £ are not related to the beam 
shape. The monopole effect arises from circular beams with unmatched main-beam full width at half maximum, 
which is described by the parameter fi = (cti — <7i)I (ci + (72), where a\ and 02 are the mean beamwidthes of first and 
second beams. The quadrupole effects arises from beams with differential ellipticities, and described by the ellipticity 
parameter e = {ox — <^y)/i'^x + fy), where ax and ay are the major and minor axes of the beam. Also, differential 
pointing, i.e. the dipole effect, is described the parameter p = pi — p2, where pi and p2 are the circular positions of 
first and second beams. These three effects can induce future spurious polarization singles from temperature leakage. 

In the formulae in Table IIIIl the functions /i , /2 and /a are experiment-specific and encapsulate the information 
about the scanning strategy which couples to the beam mismatch parameters to generate spurious polarization. The 
exact definitions of /i, /2 and /a are given in Eq. (27) in {S^, i.e. 



where 



/(m, n) ^ (e-(2rn+„)a^ , h±^^ [/(m, n) - /(to ± 2, n) (e*^'")] , D^l- (e^^) {e-^^- 



f(rn, n) and h±{'m,n) are the Fourier transformations of /(m, n) and h±{m,n). Exact definition of a can be found 
in [59| . Angular brackets represent average over measurements of a single pixel, averaged over time. In general, these 
functions are spatially-anisotropic but for simplicity, and to obtain a first-order approximation, we consider them 
constants in general. 

In the real data analysis, these generated CMB power spectra may be treated as a part of the signals, as well as 
the real signals generated by the perturbations fields, i.e. 

Cj^Cj+Cl^^,. (31) 

Thus the estimator of this contamination Dj becomes biased for the true power spectra Cj ^ by a term ^. And 
this biased estimator will induce the biased estimator for the cosmological parameters in the likelihood analysis. As 
in general, we can use the rj^^ and ntwL as the best estimator for the parameters r* and nt- Following the previous 
works [G^, we define the bias of the tensor-to-scalar ratio 5r and the spectral index 6nt as follows, 

= (^ml) - ^"^t = {ntML) - nt. (32) 

Given the beam systematics the bias of the tensor-to-scalar ratio 5r and the spectral index 5nt can be calculated by 
the following formulae (similar to the previous works [60j). 

where = Cj {sys) / a j^y . Notice that, for the requirement of the beam systematics of CMBPol mission (for instant, 
the parameter g, associated with the differential gains, satisfies g <€. 0.01% the power spectra generated by the 

beam systematics are expected to be much smaller than those generated by gravitational waves with r > 10~^ but 
the very large multipole range, where the noises are dominant. So beam systematics cannot change the value of the 
best-pivot multipole. 

Considering the CMB experiment with multi-frequency channels, the effective combined noise power spectra in Eq. 
(PS)) can be extended to the follows, considering the contribution of the systematics, 

[Nls,i]-' = E i^A') + + ^sysM ' (34) 
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where i runs though the channels. Throughout this section, we shall assume the optimistic foreground removal with 
the residual factor cr^s = q.OI for both EPIC-2m and EPIC-LC missions. We should remember that, similar to the 
previous discussion, the B-mode contamination (cr'°"^ — 0.5 for EPIC-2m, and cr''^"*' — 1 for EPIC-LC) by weak 
lensing effect is also considered throughout this section. The total CMB power spectra generated beam systematics 
can be calculated by C/(sys) = A^^^ ^(with sys) — N^ff ^(no sys). In the simplest case with single frequency channel, 
this term returns to C|'(sys) = C^^ g. 

Let us separately investigate the five systematical effects. Fig. [TO]shows C^(sys) for different values of the parameter 
g = 0.002%, 0.005% and 0.01%. In all these figures, /i = 2tt is used as the worst case. The values of i{i + l)Cf (sys) 
only weekly depends on the multipole £. As long as r > 0.01, we find C^(sys) are all smaller than those of the 
signals Cf or the noises iV^,^ ^ in the range £ < 200. In Fig. [TTl we plot the values of biases 5r and Srit induced by 
the differential gains. We find that the values of 5r and Snt strongly depends on the value of the parameter g. A 
larger g follows a larger bias. However, the uncertainties Ar and Ant are nearly independent of g unless the value 
of g is too large. We also find that, given the parameter g the value of the ratio Sr/ Ar also strongly depend on the 
tensor-to-scalar ratio r. For example, the EPIC-2m experiment with g = 0.01%, 5r is larger than Ar when r < 0.04, 
the bias is very obvious. However when r > 0.04, we have Sr < Ar, the bias is smaller than the uncertainty. 

Similar to the previous work [60| . we can define the critical value gc, which is the largest value of g as long as the 
the condition Sr/Ar < 0.1 is satisfied. In Table ITVl we list the critical values of gc for EPIC-2m and EPIC-LC, where 
r = 0.001, 0.01 and r = 0.1 are considered. Since Fig. [TT] shows that for a given g value, both Ar and 5r increase 
with the increasing of r. So the tendency of the critical gc for different input r is not trivial. After analysis, we find 
that for the assumed g value, with the increasing of r, if the increasing of Ar is more rapid than that of (5r, thus the 
case with larger r corresponds to a larger gc- This is clearly shown in Table HVl On the other hand, if the increasing 
of 5r is more rapid than that of Ar, thus the case with larger r corresponds to a smaller gc- Similar discussion is also 
applied to the nt case, as well as the cases for the other four systematics contaminations. From Table HVl we find the 
severest constraints are obtained from the requirement of the small r case. The requirement for EPIC-2m is quite 
close to that of EPIC-LC. 

Let us discuss the effect of beam gain on the determination of spectral index nt- From Eq. ([33]), we find that the 
bias ej in the lower multipole range £ < contributes a negative Snt, and the bias ej in the high multipole range 
£ > it contributes a positive 6nt. These two components are cancelled by each other and total bias Snt is expected 
to be very small, which is clearly shown in Fig. 1111 (right panel). Comparing with the bias of r, the bias Snt is much 
smaller than that of Ant. So the constraint on g obtained from the requirement of the nt is much larger than that 
from the parameter r. 

Now, let us turn to the effect of the differential monopole effect. From the formulae in Table IIIII we know the 
power spectra generated by differential monopole effect strongly depend on the beam size. Larger beam size follows 
the larger power spectra C^yg^e- This can be seen clearly in Fig. [12] Given fi = 0.1%, the value of C^(sys) is much 
larger in EPIC-LC than that in EPIC-2m. So in order to achieve a same value of Sr/ Ar = 0.1, the requirement for 
EPIC-LC is much severer than that for EPIC-2m, which are clearly shown in Table IIVI Fig. [T3] shows that, for a 
given /i, a larger r follows a larger ratio value Sr/Ar, which is correct for both EPIC-2m and EPIC-LC. In order 
to keep the tensor-to-scalar ratio in the range r g (0.001,0.1) unbiased, we should have fic — 0.029% for EPIC-2m, 
and /ic = 0.005% for EPIC-LC, where /i = 27r is adopted, we can also discuss the effect of differential monopole on 
the determination of spectral index nt- From Fig. [T^] we find that ^ is sharply peaked at the high multipole, 
where £ > £*. Thus the total contribution to Snt is always positive (see Fig. [T3l) . especially when the value of /i is 
not too small, the effect of differential monopole at the small scale is very important, and follows a fairly large bias 
for spectral index. Let us define the value of fic, where Snt/ Ant = 0.1 is satisfied. From Table HV] we find that the 
constraint on is a little severer from the parameter nt than that from the parameter r. This table shows that, the 
most severe constraint on fj, is obtained from the requirement of nt in the case of r = 0.1. 

In Figs. [T3] and [T5] we show the effect of differential pointing on the determination of gravitational waves, where 
/2 = 2tt is adopted. These figures show that similar with the case of differential monopole effect, the function C^(sys) 
generated by differential pointing is also sharply peaked at the high £, which follows that the bias of the spectral 
index Snt is positive. Comparing the left panel with the right panel in Fig. I14[ we find that for a given p, the effect of 
the differential pointing is more important for the EPIC-2m, due to the smaller instrumental noises. So more severe 
constraint on the differential pointing is followed for the EPIC-2m. In Table ITVl we find that for both EPIC-2m and 
EPIC-LC, the most severe on p comes from the requirement of parameter rif at r = 0.1. 

We have also investigated the effect of the differential quadrupole in Figs. [T^] and [T7] The effects are similar to 
those of the differential monopole. We find that EPIC-LC needs the much stricter requirement than EPIC-2m. For 
each experiment, the most severe constraint on the parameter e is obtained from the requirement of nt in the case 
with r = 0.1. 

At last, we shall discuss the effect of the differential rotation, and the results are shown in Figs. [T51 and [T51 We 
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TABLE III: The leading order contributions of tlie s;5Astematic effects to the CMB power spectra, assuming the underlying 
sky is not polarized (except for the rotation signal) 5!||, where z = (la)^% and cr is mean beamwidth of the beam, which is 
calculated by ct = SF/%/81n2. ce = cos(^), where Q is the angle between ellipse major axis of the elliptical gaussian beam 
and the horizontal i-axis of the fixed focal plane, = cos(i/)) and = sin('(/i), where ^ is the angle between the axis of 
polarization sensitivity and the major axis of the elliptical beam. A clear show of the angles Q and t^) can be found in Fig. 1 in 
[59|. The definitions of the other parameters can be found in the text, see also [s^ . 
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FIG. 10: The contribution of differential gain to the B-polarization, comparing with the signal of Cf for different r, the noise 
of CMBPol mission, and the lensed _B-polarization. Left panel is for EPIC-2m, and the right panel is for EPIC-LC. 

find that for a given e, the ratios Sr/Ar and Snt/Arit only weakly depend on the tensor-to-scalar ratio. And the 
ratio for EPIC-2m is a little smaller than that of EPIC-LC. In order to keep the parameters r and rit unbiased, the 
requirement e < 0.09° is needed for EPIC-2m, and e < 0.15° is needed for EPIC-LC. 

As a conclusion, by analyzing the effects of the systematics on the determination of gravitational waves, we find 
that the requirement of nt unbiased follows the similar or even more severe constraints for the beam systematical 
parameters. For the effects of differential monopole, pointing and quadrupole, a larger r follows a more severe 
constraint for the systematics. The critical values for the systematical parameters are listed in Table IIVI where the 
bold entries denote the most severe constraint in each case. We also find that comparing with EPIC-2m, the low cost 
EPIC-LC experiment has a much high requirement for the systematical parameters fi and e. 

VII. CONCLUSION 

The proposed CMBPol mission is the next generation of the space-based CMB experiment, which will survey the 
full sky and has the much smaller instrumental noises than Planck satellite. As one the most important tasks of this 
mission, detecting relic gravitational waves will be achieved if the tensor-to-scalar ratio r > 0.001, which will provide 
a great opportunity to study the physics in the early Universe, especially in the inflationary stage. 
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FIG. 11: The values of Sr (left panel) and 5nt (right panel) from the differential gain for different tensor-to-scalar ratio r. For 
the comparison, we plot the corresponding Ar and Ant in solid lines in the panels. 
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FIG. 12: The contribution of monopole effect to the B-polarization, comparing with the signal of Cf for different r, the noise 
of CMBPol mission, and the lensed _B-polarization. Left panel is for EPIC-2m, and the right panel is for EPIC-LC. 

In this paper, we have detailedly discussed the detection of relic gravitational waves by focusing on the constraints 
of the parameters r and the spectral index nt , which are always used to describe the primordial power spectrum of 
gravitational waves. In our discussion, we deeply investigate various contaminations for the detection, including the 
instrumental noises of CMBPol mission, the cosmic lensing contamination, the foreground contaminations, and the 
effect of various beam systematics. We found that the cosmic lensing becomes the dominant noise sources, comparing 
the instrumental noises for both EPIC-2m and EPIC-LC projects. Different from Planck satellite, if r > 0.01, the 
detection of gravitational waves mostly depends on the observation at multipole £ ~ 100, the peak of the i3-mode 
polarization. However, the reionization peak at £ ~ 10 still plays a crucial role for the determination of spectral index 
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FIG. 13: The values of Sr (left panel) and 5nt (right panel) from the monopole effect for different tensor-to-scalar ratio r. For 
the comparison, we plot the corresponding Ar and Ant in solid lines in the panels. 




FIG. 14: The contribution of differential pointing to the B-polarization, comparing with the signal of Cf for different r, the 
noise of CMBPol mission, and the lensed _B-polarization. Left panel is for EPlC-2m, and the right panel is for EPIC-LC. 



rif. We also found that if the foreground contaminations cannot be well controlled, the reionization peak may be 
unobservable, which could deeply increase the uncertainty of rit- At the same time, we have investigated the effect 
of various beam systematics on the detection of gravitational waves, which mainly cause a bias on the cosmological 
parameters. In order to keep these biases small enough, the requirements for the beam systematical parameters are 
quite severe, especially for the EPIC-LC mission. 
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FIG. 15: The values of Sr (left panel) and Snt (right panel) from the differential pointing for different tensor-to-scalar ratio r. 
For the comparison, we plot the corresponding Ar and Ant in solid lines in the panels. 
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FIG. 16: The contribution of quadrupole effect to the B-polarization, comparing with the signal of Cf for different r, the noise 
of CMBPol mission, and the lensed B-polarization. Left panel is for EPIC-2m, and the right panel is for EPIC-LC. 
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FIG. 17: The values ol Sr (left panel) and 5nt (right panel) from the quadrupole effect for different tensor-to-scalar ratio r. 
For the comparison, we plot the corresponding Ar and Ant in solid lines in the panels. 




FIG. 18: The contribution of differential rotation to the _B-polarization, comparing with the signal of Cf for different r, the 
noise of CMBPol mission, and the lensed _B-polarization. Left panel is for EPIC-2m, and the right panel is for EPIC-LC. 



[1] J. Bock et al, Task force on cosmic microwave background research, arXiv:astro-ph/0604101 . 

[2] L. P. Grishchuk, Zh. Eksp. Teor. Fiz. AmpUfication of gravitational waves in an isotroptc universe, 67, 825 (1974) [Sov. 
Phys. JETP 40, 409 (1975)]; Graviton creation in the early Universe, Ann. NY Acad. Sci. 302, 439 (1977); Primordial 



22 




FIG. 19: The values of Sr (left panel) and Srit (right panel) from the differential rotation for different tensor-to-scalar ratio r. 
For the comparison, we plot the corresponding Ar and Ant in solid lines in the panels. 



TABLE IV: Systematics tolerance for EPIC- 2m and EPIC-LC, where different nominal values of r are considered. In each box, 
the left value is for the parameter r and the right value is for the parameter rit. 





Nominal value 


1% V 2ir 


1% V 2ir 


1" V 2vr 


1% 


ec[deg] 


EPIC-2m 


r=0.001 
r=0.01 
r=0.1 


0.0017 & 

0.0021 & 

0.0031 & 


0.061 & 0.042 
0.049 & 0.030 
0.029 & 0.019 


0.096 & 0.085 
0.087 & 0.071 
0.068 & 0.054 


0.87 & 0.59 
0.70 & 0.43 
0.41 & 0.27 


0.11 & 0.21 
0.11 & 0.14 
0.09 & 0.10 


EPIC-LC 


r=0.001 
r=0.01 
r=0.1 


0.0020 & 

0.0026 & 

0.0039 & 


0.0073 & 0.0061 
0.0064 & 0.0053 
0.0050 & 0.0040 


0.15 & 0.14 
0.14 & 0.13 
0.13 & 0.11 


0.10 & 0.09 
0.09 & 0.08 
0.07 & 0.06 


0.15 & 0.22 
0.17 & 0.19 
0.18 & 0.16 



gravitons and possibility of their observation, Pis'ma Zh. Eksp. Teor. Fiz. 23, 326 (1976) [JETP Lett. 23, 293 (1976)]. 
[3] L. P. Grishchuk, Discovering relic gravitational waves in cosmic microwave background radiation, Chapter in the "General 
Relativity and Johh Archibald Wheeler", edited by I. Ciufolini and R. Matzner, (Springer, New York, 2010, pp. 151-199), 
arXiv:0707.3319 

[4] A. A. Starobinsky, Spectrum of relict gravitational radiation and the early state of the Universe, JETP Lett. 30, 682 
(1979). 

[5] V. A. Rubakov, M. Sazhin and A. Veryaskin, Graviton creation in the inflationary universe and the grand unification 

scale, Phys. Lett. B 115, 189 (1982). 
[6] S. Sasaki, Large scale quantum fluctuations in the inflationary universe, Prog. Theor. Phys. 76, 1036 (1986). 
[7] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Theory of cosmological perturbations, Phys. Rep. 215, 203 

(1992). 

[8] L. P. Grishchuk, Relic gravitational waves and their detection. Led. Notes Phys. 562, 167 (2001); 

Y. Zhang, Y. F. Yuan, W. Zhao and Y. T. Chen, Rehc gravitational waves m the accelerating Universe, Class. Quant. 
Grav. 22, 1383 (2005); 

T. L. Smith, H. V. Peiris and A. Cooray, Deciphering inflation with gravitational waves: Cosmic microwave background 
polarization vs direct detection with laser interferometers, Phys. Rev. D 73, 123503 (2006); 
W. Zhao and Y. Zhang, Relic gravitational waves and their detection, Phys. Rev. D 74, 043503 (2006). 
[9] Y. Watanabe and E. Komatsu, Improved calculation of the primordial gravitational wave spectrum in the standard 
model, Phys. Rev. D 73, 123515 (2006); 

L. A. Boyle and P. J. Steinhardt, Probing the early universe with inflationary gravitational waves, Phys. Rev. D 77, 063504 
(2008); 

M. Giovannini, Production and detection of relic gravitons in quintessential inflationary models, Phys. Rev. D 60, 123511 



23 



(1999), Spikes in the relic graviton background from quintessential inflation, Class. Quant. Grav. 16, 2905 (1999), Stochastic 
backgrounds of relic gravitons, TACDM paradigm and the stiff ages, Phys. Lett. B 668, 44 (2008), Thermal history of the 
plasma and high-frequency gravitons, Class. Quant. Grav. 26, 045004 (2009). 
[10] D. H. Lyth and A. Riotto, Particle physics models of inflation and the cosmological density perturbation, Phys. Rep. 314, 
1 (1999). 

[11] S. Kachru et ai. Towards inflation in string theory, JCAP, 0310, 013 (2003); 

D. Baumann and L. McAllister, A microscopic limit on gravitational waves from D-brane inflation, Phys. Rev. D 75, 
123508 (2007); 

E. Silverstein and D. Tong, Scalar speed limits and cosmology: Acceleration from D-ceeleration, Phys. Rev. D 70, 103505 
(2004); 

R. Kallosh and A. Liddle, Testing string theory with CMB, JCAP, 0704, 017 (2007). 
[12] A. Polnarev, Polarization and anisotropy induced in the microwave background by cosmological gravitational waves, Sov. 
Astron. 29, 607 (1985); 

D. Harari and M. Zaldarriaga, Polarization of the microwave background in inflationary cosmology, Phys. Lett. B 319, 96 
(1993); 

L. P. Grishchuk, Cosmological perturbations of quantum-mechanical origin and anisotropy of the microwave back- 
ground, Phys. Rev. Lett. 70, 2371 (1993); 

R. Crittenden, J. R. Bond, R. L. Davis, G. Efstathiou and P. J. Steinhardt, Imprint of gravitational waves on the cosmic 
microwave background, Phys. Rev. Lett. 71, 324 (1993); 

R. A. Frewin, A. G. Polnarev and P. Coles, Gravitational waves and the polarization of the cosmic microwave back- 
ground, Mon. Not. R. Astron. Soc. 266, L21 (1994). 
[13] U. Seljak and M. Zaldarriaga, Signature of gravity waves in the polarization of the microwave background, Phys. Rev. Lett. 
78, 2054 (1997); 

M. Kamionkowski, A. Kosowsky and A. Stebbins, A probe of primordial gravity waves and vorticity, Phys. Rev. Lett. 78, 
2058 (1997). 

[14] M. Zaldarriaga and U. Seljak, All-sky analysis of polarization in the microwave background, Phys. Rev. D 55, 1830 (1997). 
[15] J. R. Pritchard and M. Kamionkowski, Cosmic microwave background fluctuations from gravitational waves: An analytic 
approach, Ann. Phys. (N.Y.) 318, 2 (2005); 

W. Zhao and Y. Zhang, Analytic approach to the CMB polarization generated by relic gravitational waves, Phys. Rev. D 
74, 083006 (2006); 

D. Baskaran, L. P. Grishchuk and A. G. Polnarev, Imprints of relic gravitational waves in cosmic microwave background 
radiation, Phys. Rev. D 74, 083008 (2006); 

T. Y. Xia and Y. Zhang, Analytic spectra of CMB anisotropics and polarization generated by relic gravitational waves with 
modification due to neutrino free- streaming, Phys. Rev. D 78, 123005 (2008), Approximate analytic spectra of reionized 
CMB anisotropics and polarization generated by relic gravitational waves, Phys. Rev. D 79, 083002 (2009). 
[16] B. G. Keating, A. G. Polnarev, N. J. Miller and D. Baskaran, The polarization of the cosmic microwave background due 
to primordial gravitational waves. Int. J. Mod. Phys. A 21, 2459 (2006); 

R. Flauger and S. Weinberg, Tensor microwave background fluctuations for large multipole order, Phys. Rev. D 75, 123505 
(2007); 

Y. Zhang, W. Zhao, X. Z. Er, H. X. Miao and T. Y. Xia, Relic gravitational waves and CMB polarization in the accelerating 
Universe, Int. J. Mod. Phys. D 17, 1105 (2008). 
[17] A. G. Polnarev, N. J. Miller and B. G. Keating, CMB temperature polarization correlation and primordial gravitational 
waves , Mon. Not. R. Astron. Soc. 386, 1053 (2008); 

N. J. Miller and B. G. Keating and A. G. Polnarev, CMB temperature polarization correlation and primordial gravitational 
waves II: Wiener filtering and tests based on Monte Carlo simulations, arXiv:0710.3651 

[18] W. Zhao, D. Baskaran and L. P. Grishchuk, On the road to discovery of relic gravitational waves: The TE and BB 
correlations in the cosmic microwave background radiation, Phys. Rev. D 79, 023002 (2009). 

[19] E. Komastu et al, Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Cosmological interpreta- 
tion, arXiv:1001.4538 

[20] W. Zhao, D. Baskaran and L. P. Grishchuk, On the road to discovery of relic gravitational waves: The TE and BB 
correlations in the cosmic microwave background radiation, Phys. Rev. D 79, 023002 (2009), Stable indications of relic 
gravitational waves in Wilkinson Microwave Anisotropy Probe data and forecasts for the Planck mission, Phys. Rev. D 80, 
083005 (2009), Relic gravitational waves in light of the 7-year Wilkinson Microwave Anisotropy Probe data and improved 
prospects for the Planck mission, Phys. Rev. D 82, 043003 (2010); 

W. Zhao and L. P. Grishchuk, Relic gravitational waves: Latest revisions and preparations for new data, Phys. Rev. D 82, 
123008 (2010). 

[21] Planck Collaboration, The science programme of Planck, ^arXiviastro-ph / 0604069] . 

[22] W. Zhao, Detecting relic gravitational waves in the CMB: Comparison of different methods, Phys. Rev. D 79, 063003 
(2009). 

[23] Y. Z. Ma, W. Zhao and M. L. Brown, Constraints on standard and non-standard early Universe models from CMB B-mode 

polarization, JCAP, 1010, 007 (2010). 
[24] W. Zhao and D. Baskaran, Detecting relic gravitational waves in the CMB: Optimal parameters and their constraints, Phys. 

Rev. D 79, 083003 (2009). 

[25] B. G. Keating et al, BICEP: A large angular scale CMB polarimeter, in Polarimetry in Astronomy, edited by Silvano 



24 



Fineschi, Proceedings of the SPIE, 4843 (2003); 

C. Pryke et al, QUaD Collaboration, Second and third season QUaD CMB temperature and polarization power spectra, As- 
trophys. J. 692, 1247 (2009); 

A. C. Taylor, Clover Collaboration, Clover - A B-mode polarization experiment, New Astron. Rev. 50, 993 (2006); 
J. E rrard, The new generation CMB B-mode polarization experiment: POLARBEAR, arXiv:1011.0763 
http:/ /bolo. berkeley.edu/polarbear/index. html 

D. Samtleben, Measuring the cosmic microwave background radiation (CMBR) polarization with O j7/ii/r. larXiv:0802.2657l 
[http:/ /quiet. uchicago.edu/ 

J.A. Rubino-Martin et al.. The Quijote CMB Experiment, iarXiv:0810.3141l 



http:/ /www. iac.es/project 


/cmb/quijote/ 


http:/ /pole. uchicago.edu/ 




http : / / wwwphy.princeton.edu / act / 1 



[26] http:/ /groups. physics. umn.edu/cosmology/ebex/index. html 

A. Kogut, et al., PAPPA: Primordial Anisotropy Polarization Pathfinder Array, New Astron. Rev. 50, 1009 (2006); 

B. P. Crill, et al., SPIDER: A balloon-borne large-scale CMB polarimeter . larXiv:080 7.1548', 

[27] W. Zhao and W. Zhang, Detecting relic gravitational waves m the CMB: Comparison of Planck and ground-based experi- 
ments, Phys. Lett. B 677, 16 (2009). 

[28] D. Baumann et al, CMBPol mission concept study: Probing inflation with CMB polarization, larXiv:0811.3919| 

[29] B-Pol Collaboration, B-Pol: Detecting primordial gravitational waves generated during inflation, Exper. Astron. 23, 5 

(2009); 

http://www.b-pol.org/index.php l 



[30] 'http: / / cmbpol.kek.jp/litebird / 
[31] L. P. Grishchuk and M. Solokhin, Spectra of relic gravitons and the early history of the Hubble parameter, Phys. Rev. D 
43, 2566 (1991); 

A. Kosowsky and M. S. Turner, CBR anisotropy and the running of the scalar spectral index, Phys. Rev. D 52, 1739 (1995). 

[32] E. D. Stewart and D. H. Lyth, A more accurate analytic calculation of the spectrum of cosmological perturbations produced 
during inflation, Phys. Lett. B 302, 171 (1993). 

[33] H. V. Peiris et al, First Year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Implications for inflation, As- 
trophys. J. Suppl. Ser. 148, 213 (2003). 

[34] E. Komatsu et al., Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Cosmological interpreta- 
tion, Astrophys. J. Suppl. Ser. 180, 330 (2009). 

[35] A. Lewis, A. Challinor and A. Lasenby, Efficient computation of CMB anisotropics m closed FRW models, Astrophys. J. 
538, 473 (2000); 
http:/ /camb.info/ 

[36] M. J. Mortonson and W. Hu, Model-independent constraints on reionization from large-scale CMB polarization, Astrophys. 
J. 672, 737 (2008); 

M. J. Mortonson and W. Hu, Impact of reionization on CMB polarization tests of slow-roll inflation, Phys. Rev. D 77, 
043506 (2008); 

M. J. Mortonson and W. Hu, Reionization constraints from five-year WMAP data, Astrophys. J. 686, L53 (2008). 
[37] S. Hamimeche and A. Lewis, Likelihood analysis of CMB temperature and polarization power spectra, Phys. Rev. D 77, 
103013 (2008). 

[38] W. J. Percival and M. L. Brown, Likelihood methods for the combined analysis of CMB temperature and polarisation power 
spectra, Mon. Not. R. Astron. Soc. 372, 1104 (2006); 

H. K. Eriksen and I. K. Wehus, Marginal distributions for cosmic variance limited CMB polarization data, Astrophys. J. 
Suppl. Ser. 180, 30 (2009). 

[39] A. Lewis and S. L. Bridle, Cosmological parameters from CMB and other data: A Monte Carlo approach, Phys. Rev. D 
66, 103511 (2002). 

[40] A. Lewis, Harmonic E/B decomposition for CMB polarization maps, Phys. Rev. D 68, 083509 (2003); 

A. Lewis, A. Challinor and N. Turok, Analysis of CMB polarization on an incomplete sky, Phys. Rev. D 65, 023505 (2001); 
E. F. Bunn, M. Zaldarriaga, M. Tegmark and A. de Oliveira-Costa, E/B decomposition of finite pixelized CMB maps, Phys. 
Rev. D 67, 023501 (2003); 

E. F. Bunn, E/B mode mmno, farXiv: 08 1 1 . 1 1 1 1 Efficient decomposition of cosmic microwave background polarization maps 
into pure E, pure B, and ambiguous components, arXiv: 1008. 0827' 
[41] K. M. Smith, Pseudo-Cl estimators which do not mix E and B modes , Phys. Rev. D 74, 083002 (2006); 

J. Grain, M. Tristram and R. Stompor, Polarized CMB power spectrum estimation using the pure pseudo-cross-spectrum 
approach, Phys. Rev. D 79, 123515 (2009); 

K. M. Smith and M. Zaldarriaga, General solution to the E-B mixing problem, Phys. Rev. D 76, 043001 (2007). 

[42] W. Zhao and D. Baskaran, Separating E and B types of polarization on an incomplete sky, Phys. Rev. D 82, 023001 (2010); 
J. Kim and P. Naselsky, CMB E/B decomposition of incomplete sky: a pixel space approach. Astronomy and Astrophysics 

519, A104 (2010); 

J. Kim, How to make clean separation of CMB E and B mode from incomplete sky data, larXiv:1010.2636) 

[43] J. Bock et al.. Study of the Experimental probe of Inflationary Cosmology (EPIC)-Intemediate mission for NASA's Einstein 
Inflation Probe, arXiy:0906.1188. 



25 



[44] M. Zaldarriaga and U. Seljak, Gravitational lensing effect on cosmic microwave background polarization, Phys. Rev. D 58, 
023003 (1998). 

[45] W. Hu and T. Okamoto, Mass reconstruction with CMB polarization, Astrophys. J. 574, 566 (2002); 

T. Okamoto and W. Hu, CMB lensmg reconstruction on the full sky, Phys. Rev. D 67, 083002 (2003). 
[46] C. M. Hirata and U. Seljak, Reconstruction of lensmg from the cosmic microwave background polarization, Phys. Rev. D 

68, 083002 (2003); 

U. Seljak and C. M. Hirata, Gravitational lensing as a contaminant of the gravity wave signal m the CMB, Phys. Rev. D 

69, 043005 (2004). 

[47] A. Lewis and A. Chalhnor, Weak gravitational lensing of the CMB, Phys. Rep. 429, 1 (2006). 
[48] K. M. Smith et at, CMBPol mission concept study: Gravitational lensing, |arXiv:0 811.3916 

[49] L. Knox and Y. S. Song, Limit on the detectability of the energy scale of inflation, Phys. Rev. Lett. 89, 011303 (2002); 
M. Kesden, A. Cooray and M. Kamionkowski, Separation of gravitational-wave and cosmic-shear contributions to cosmic 
microwave background polarization, Phys. Rev. Lett. 89, 011304 (2002). 

[50] C. L. Bennett et al. First Year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Foreground emission, As- 
trophys. J. Suppl. Ser. 148, 97 (2003). 

[51] J. Dunkely et al, CMBPol mission concept study: Prospects for polarized foreground removal, larXiv:0811.3915] 

[52] L. Verde, H. Peiris and R. Jimenez, Considerations m optimizing CMB polarization experiments to constrain inflationary 
physics, JCAP, 0601, 019 (2006). 

[53] L. Page et al.. Three Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Polarization Analysis, Astro- 
phys. J. Suppl. Ser. 170, 335 (2007). 

[54] E. M. Finkbeiner et al.. Extrapolation of galactic dust emission at 100 microns to CMBR frequencies using FIRAS, 
Astrophys. J. 624, 10 (2001). 

[55] J. Kim, P. Naselsky and P. R. Christensen, CMB polarization map derived from the WMAP 5 year data through the 
harmonic internal linear combination method, Phys. Rev. D 79, 023003 (2009); 

G. Efstathiou, S. Gratton and F. Paci, Impact of galactic polarized emission on B-mode detection at low multipoles, Mon. 
Not. Roy. Astron. Soc. 397, 1355 (2009). 



[56] W. Zhao and Q. G. Huang, Testing inflationary consistency relations by the potential CMB observations, larXiv: 1 101 . 31631 

[57] W. Hu, M. M. Hedman and M. Zaldarriaga, Benchmark parameters for CMB polarization experiments, Phys. Rev. D 67, 
043004 (2003). 

[58] M. Su, A. P. S. Yadav and M. Zaldarriaga, Impact of instrumental systematic contamination on the lensing mass recon- 
struction using the CMB polarization, Phys. Rev. D 79, 123002 (2009); 

M. Su, A. P. S. Yadav, M. Shimon and B. G. Keating, Impact of instrumental systematics on the CMB bispec- 
iraw, i,arXiv:10 10.1957 

[59] M. Shimon, B. Keating, N. Ponthieu and E. Hivon, CMB polarization systematics due to beam asymmetry: Impact on 
inflationary science, Phys. Rev. D 77, 083003 (2008); 

N. J. Miller, M. Shimon and B. G. Keating, CMB beam systematics: Impact on lensing parameter estimation, Phys. Rev. 
D 79, 063008 (2009); CMB polarization systematics due to beam asymmetry: Impact on cosmological birefringence, Phys. 
Rev. D 79, 103002 (2009). 

[60] D. O'Dea, A. Challinor and B. R. Johnson, Systematic errors in cosmic microwave background polarization measure- 
ments, Mon. Not. Roy. Astron. Soc. 376, 1767 (2007). 

[61] Note that in the real analysis throughout this paper, similar to [2^, we have not considered the frequency channels at 
30GHz and 340GHz proposed in 28] for EPIC-2m, and the frequency channels at 30GHz and 300GHz for EPIC-LC. 

[62] Note that, in Sec.|Vl we will show that the total noise level will be increased if the reduced foreground contaminations are 
considered. However, the noise level with high multipole £ > 100 nearly keeps same. 




26 



